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Motivated by recent tunnelling and Andreev-reflection experiments, we study the conditions for 
a quantum transition within the superconducting phase of the cuprates, in which a bulk imaginary 
(time-reversal breaking) id^y component appears in addition to the d^2_yi order parameter. 

We examine in detail the role of some important physical features of the cuprates. In particular we 
show that a closed Fermi surface, a bilayer splitting, an orthorhombic distortion, and the proximity 
to a quantum critical point around optimal doping favor the appearance of the imaginary component. 
These findings could explain why the mixed d^2_y2 + id^y order parameter is observed in YBCO 
and LSCO, and suggest that it could appear also in Bi2212. We also predict that, in all cuprates, 
the mixed state should be stable only in a limited doping region all contained beneath the d^2_y2 
dome. The behavior of the specific heat at the secondary transition is discussed. 

PACS numbers: 74.30.Rp, 74.20.Fg, 74.25.Dw 



I. INTRODUCTION 

The evidence for a d^2_y2 symmetry of the supercon- 
ducting pairing in the High-Tc Cuprate Superconductors 
(HTCS) is nowadays overwhelmingi. Nonetheless it is 
well known that a secondary gap component of different 
symmetry can develop either spontaneously or driven by 
external factors, like a magnetic field, magnetic impuri- 
ties as well as the proximity to an interfac o^i'^i'^i^ . The 
first important consequence is that the nodal lines of the 
gap located along the diagonals of the Brillouin zone (FX 
direction) can be removed, deeply altering the low energy 
properties. Secondly, this kind of superconducting state 
can break the time-reversal symmetry with important 
consequences on the physical properties^. Specifically a 
superconducting state with mixed dx2-y2 + idxy symme- 
try has been invoked in Refs.^'^'^ to explain the thermal 
conductivity anomalies in Bi2Sr2CaCu208 (Bi2212)^'^, 
in RefiiS to account for the properties of single quasipar- 
ticle states observed in the vortex core^^, and in Ref.^^ 
as a possible explanation of the momentum and tempera- 
ture dependence of the damping of quasiparticles. Scan- 
ning tunnelling experiments in Yi_j,Caj,Ba2Cu307_:r 
(YBC0)iii4 and in La2-:.Sr^Cu04 (LSCO)i^ have also 
shown clear deviations from the pure (ia;2_j,2-wave sym- 
metry. While the LSCO results have been originally in- 
terpreted in terms of a d -I- s wave, a recent reanalysis 
based on the results of the present report has shown that 
the d,j.2_y2 + idxy symmetry is a more faithful description 
of the experiments^. 

Our main motivation comes indeed from tunnelling ex- 
periment on YBCO films^'*, which have evidenced that 
the zero-bias conductance peak associated with a pure 
dx^~y^ gap function^'' undergoes a doping- and magnetic- 
field- dependent splitting. This can be interpreted^* in 
terms of a secondary time-reversal breaking component 
of the superconducting gap. It is quite important to no- 



tice that the amplitude of the splitting has a power-law 
dependence on doping and magnetic-field. The (zero- 
field) splitting is in fact linear in the deviation from the 
optimal doping for overdoped samples, and it is zero for 
underdoped samples. Moreover, starting from a finite 
(zero) value at zero field, the splitting increases linearly 
in the overdoped (underdoped) materials with the mag- 
netic field, while it displays a square-root dependence 
close to optimal doping. These power-law dependen- 
cies are strong indications of a critical behavior due to a 
second-order phase transition and are hardly accounted 
for by surface or disorder effects alone. The above result 
suggests therefore a doping-driven bulk transition. The 
observed fraction of d^y pairing increases with S, but it 
must be noticed that it never exceeds the 20% of the 
principal one in all the measured samples^. These ex- 
perimental situation leaves the door open for at least two 
scenarios. In the first one, discussed in a different context 
in Ref.^§i, the dxy component keeps increasing by further 
doping and eventually the gap becomes of pure dxy sym- 
metry. In the alternative scenario, which is supported by 
our analysis, the dxy fraction does not monotonically in- 
crease, but it starts to decrease at some doping, leading 
to a smaller dome beneath the dx2-y2, where the mixed 
order parameter is stable. Our scenario is also consistent 
with other experiments which put severe upper bounds 
for any imaginary component of the gap (see Ref.-- for a 
recent discussion). 

The experimental framework described above sug- 
gested us to carry out a systematic analysis of the specific 
properties which may favor the occurrence of a mixed su- 
perconducting state in some cuprates and not in others. 
To achieve this goal, we quantitatively analyze the role of 
band-structure effects (Van Hove singularity (VHS), bi- 
layer splitting, and orthorhombic distortion), and of the 
form of the pairing interaction in determining the weight 
of the secondary component (if any), the features of the 
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phase diagram and the behavior of observables Hke the 
specific heat. We rely on the simple BCS approach to 
make the role of the above ingredients more transparent. 
We notice that the fact that the d^2_y2 + id^y order pa- 
rameter, when observed'^^, establishes in the overdoped 
regime, where strong correlation effects are less impor- 
tant. 

In Sec. II we introduce the formalism and discuss in 
different subsections the effects of the Fermi surface (A), 
the bilayer splitting (B), the orthorhombic distortion (C) 
and of a QCP (D). In Sec. Ill we discuss how a phase 
diagram for the secondary component can be drawn, and 
Sec. IV is devoted to the conclusions. 



II. METHOD AND RESULTS 

General theoretical constraints for the appearance of a 
bulk mixed order parameter have already been studied, 
e.g., in Refsi^iSS. In particular, in the first of these works 
it has been shown that the mixed state is more likely to 
establish close to a surface than in the bulk and that the 
idxy wave is the most favored secondary component in 
both cases. 

Here we recall some basic symmetry arguments and 
derive the BCS equations. The strong anisotropy of the 
HTCS allows us to focus on a purely two-dimensional lat- 
tice system, characterized by the Civ point group. Only 
the four harmonics corresponding to the one-dimensional 
irreducible representations of the group are compatible 
with singlet pairing. We can write them as^ 
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Then, the gap function can be written as 



(1) 



where the index 77 = s+ , s ,d'^,d labels the different 
representations and r = {x,y), with integer x,y, denotes 



a lattice site (we take the lattice spacing 



1). 



The primed sum is restricted to inequivalent sites under 
the symmetry of the lattice (i.e., to the different lattice 
distances). If we require the invariance of |Akp under 
the point group transformations, the gap function has 
to transform either like a single representation, or like 
a complex (time-reversal breaking) combination of the 
form + jA^ with 77 7^ C- 

Expanding the BCS equation in the same harmonics, 
it is easy to realize that the appearance of a harmonic 
labeled by r is controlled by the correspondent compo- 
nent of the pair potential V{r). In a generic model for 



the cuprates the on-site interaction (Vq) is repulsive (pos- 
itive) due to strong local repulsion. The longer distance 
part (Vi^o) is instead attractive (negative) and it is rea- 
sonably assumed to be a decreasing function of distance. 
We can therefore restrict to the smallest distance har- 
monics. In the following we always use the symbol Vi 
to denote the modulus of the V{r = i). For r = 0, the 
only nonzero harmonic belongs to the s"*" representation, 
and corresponds to the isotropic s-wave, which is ruled 
out by the repulsive Vq. For r = (±1, 0); (0, ±1) = 1 
the s"*" and d'^ representations are associated with the 
extended 5^:2+1,2 and dx2_y2 waves, respectively. For 
r = (±1,±1) = 2, the non- vanishing harmonics are re- 
lated to the Sxy and d^y waves (s+ and d^ representa- 
tions). 

If we take a complex combination of the form Aiw^.^^ + 

iA2Z«2.k: with both Ai and A2 real, each gap parameter 
affects the other only through the quasiparticle spectrum 
and the secondary pairing can be viewed as a BCS cou- 
pling between the Bogoljubov quasiparticles of the pri- 
mary superconducting phase. The BCS equations then 
read 
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where (3 is the inverse temperature, = [^^ -f 

[A^wl^f -f {A^wi.^fY'^ with ry = C = 5+,^" 

(77 ^ C). The band dispersion is ^k = —2t{coskx + 
cosfcy) -I- At' coskx cos ky — fi where t and t' are nearest 
and next-nearest neighbor hoppings and fj, is the chem- 
ical potential. Density functional theory gives a typical 
value for the hopping t in the HTCS of 200 meV, while 
t' assumes different values in the various compounds^. 
These values can however be reduced by a factor of two 
by correlation effects, that can be estimated by fitting 
the angular resolved photoemission spectra with a tight- 
binding model. As we will see later in more detail, the 
value of t' has important effect on the DOS, since it con- 
trols the position of the Van Hove singularity (VHS). 

Being the largest attractive coupling Vi controls the 
principal component of the gap, which can be dx2-y2 
or Sx2j^.y2. For typical values of the particle density 
in the HTCS, both low doping and the VHS of the 
two-dimensional DOS strongly favor the dx2-y2 sym- 
metry, which thus represents the principal component 
of the superconducting gap even in the simplified BCS 
approach^^. We take henceforth 77 = d+ in Eqs. (O,®- 
The spectrum for the secondary <^ component is therefore 
(pseudo)gapped by the principal dx2-y2 gap Ai. As a 
consequence, the mixed order parameter establishes only 
for V2 larger than a critical value ¥2"^ , as opposed to the 
case of the Cooper instability in the metallic phase, which 
takes place for arbitrarily small coupling. The critical 
coupling for both dxy and Sxy secondary pairing is deter- 
mined by solving the BCS equations Q , Q with A2 = 
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and C = d , s+ respectively. The only difference between 
these two cases is the different form of the harmonic W21C 
that weights the same kernel in 12- This kernel, in turn, is 
affected by the presence of the d^2_y2 wave. The dx2_y2 
gap in the i?k spectrum strongly suppresses the contri- 
butions from the M points (0, ±7r), (±7r, 0), while it does 
not affect much the nodal lines. The d^y wave has nodes 
along the FM directions and is maximum along the diag- 
onals, while the Sxy has its maximum contributions from 
the M points, just like the dx2-y2. It is therefore clear 
that the dxy is most likely candidate for the secondary 
component^. 

A numerical solution of the above Eqs. (0) , Q confirms 
the above arguments, finding always larger for s^y 
than for d^y Moreover, the appearance of the s^y com- 
ponent is generally associated to a first-order transition 
from pure dx2-y2 to pure s^y, since the continuous mixing 
is unlikely due to the strong competition between the two 
harmonics. This kind of scenario is clearly incompatible 
with the experiments, in which the secondary component 
appears in a continuous way. The continuous transition 
from dx^-y^ to d^^-y^ + idxy gap is therefore the most 
natural candidate if a pure _y2 is to be modified by a 
small secondary component, as suggested by the experi- 
ments on YBCOi'^ and on LSCOi^. 

Next, we evaluate the critical value of V2 to get dxy 
pairing, V^^^ = V^"^ , assuming a given value of Vi. A 

similar analysis has been performed in Ref^ for a contin- 
uous model, and in Ref.^° for a lattice model. In the lat- 
ter case, in which the effect of a two-dimensional DOS is 
considered, it is always found that > V\. This would 
imply that a secondary component is possible only tak- 
ing a potential which increases with the distance, at least 
going from nearest to next-nearest neighbors. This con- 
dition seems hard to be realized in a microscopic model, 
where the interaction is naturally a decreasing function of 
distance. Similar results have been obtained also within 
a t-J\-J2 model iniSi, where the dx2-y2 -I- idxy and pure 
dxy estabhsh, respectively, for J2I J\ — 1-4 and 2.2, which 
are clearly not representative of the cuprates. In the case 
of RefiiSij V'^'^ /Vi may be smaller than one for some values 
of the parameters, but this only comes from the assump- 
tion of a flat DOS, which is clearly not representative of 
the HTCS. 

Having in mind that for a plausible microscopic poten- 
tial V2 < Vi, in the following we explore the conditions for 
the appearance of the secondary harmonic obeying this 
constraint. First we notice that the larger is the value of 
Vi is, the smaller is the ratio V'^'^/Vi. In other words, if 
both the couplings are large, a secondary component may 
appear also for V2 <V\- The above claim can be under- 
stood as follows: If we increase Vi, Ai also increases, and 
it has different effects on the two equations. The integral 
I2 in Eq. does not change much by increasing Ai, 
due to the w'^.^^ factor, while 1\ instead contains w^^, 
and it is therefore substantially reduced by Ai. Since, 
from Eqs. V'' jVx = h/h, this implies that, at 

fixed band dispersion, the larger is Vi, the smaller is the 



ratio V'^^'/Vi and then the mixed state is more favored. 

The argument above shows that strong coupling is nec- 
essary for the appearance of a secondary componenti24. 
Nevertheless, Vi can not be increased indefinitely with- 
out pushing the value of Ai to unphysical values. If we 
want to increase the (i^2_j^2 -coupling at fixed Ai we can 
introduce a cut-off loq in the integrals. This cut-off is 
not just an artifact to push the system to stronger cou- 
pling, but naturally measures the typical energy scale of 
the attractive potential. In the forthcoming calculations 
we take ujq — 50 meV = 0.25tm^. As we see in the fol- 
lowing, this choice of cut-off leads us to reasonable phase 
diagrams without assuming unphysical values of the gap. 

Now we discuss the role of different physical effects on 
the value of the critical V2 ■ Some preliminary consider- 
ations may be done before numerically solving the BCS 
equations. In a BCS approach, the bandstructure influ- 
ences the values of the gap and of the critical tempera- 
ture only through the DOS. In particular, the VHS plays 
the main role, strongly favoring the pure dx2-y2 pair- 
ing. Other harmonics may be stabilized only for doping 
values for which the chemical potential is sufficiently far 
from the VHS. This means that, if we want the dxy to 
appear for dopings, like the ones of the experiments, the 
VHS must be relatively close to half-filling. This is re- 
alized for cuprates with closed Fermi surface like LSCO. 
To explain the appearance of the secondary component 
in YBCO we therefore have to invoke some other effects, 
like orthorhombic distortion or bilayer splitting. Both 
this effects are expected to help the secondary compo- 
nent to appear, because they "regularize" the VHS. We 
discuss the role of these effects in the following subsec- 
tions. The final subsection is instead dedicated to the 
form of the pairing potential V{r), and in particular to 
the possible role of a QCP in the charge sector close to 
optimal dopingSS.. 



A. Shape of the Fermi Surface 

As we mentioned above, the position of the VHS plays 
an important role in determining the symmetry of the or- 
der parameter within the BCS approach. The ratio t' /t 
is a direct measure of the position of the VHS, and conse- 
quently of the shape of the Fermi surface. For t' — 0, the 
VHS lies at the chemical potential at half-filling. Increas- 
ing t' /t, the density for which the chemical potential lies 
at the VHS moves away from half-filling. Consequently, 
the stability of the pure dx2-y2 wave extends to larger 
doping, making the appearance of the dxy more difficult. 
We numerically solve Eqs. Q and © to make the anal- 
ysis more quantitative. The results are reported in Fig. 
1. Here we take Vi constant as a function of doping, just 
to avoid the complications of many varying parameters 
and to disentangle the effect of the VHS. It is understood 
that Vi will depend on doping in a realistic description 
of the cuprates. 

According to the discussion above (V2 < Vi) , the sec- 
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FIG. 1: V'' jVx as a function of the doping 5 for Vxjt = 
1.15, ujQjt — 0.25 and different values of the ratio t' /t. When 
the curve crosses the line V^^ /Vi = 1 the appearance of the 
secondary harmonics becomes possible with V2 <Vi. 



FIG. 2: V'"'/Vi as a function of the doping S for Vi/t = 1.15, 
ujo/t — 0.25 and different values of the splitting parameter 
t± and of the orthorhombic parameter y = (ty ~ tx)/tx- The 
upper curve is the same as that shown in Fig0 



ondary harmonics can only appear, roughly speaking, 
when V^^ /Vi < 1. The appearance of the secondary 
harmonic is more likely as the doping increases, but the 
minimal doping 6min above which V^^ /Vi < 1 increases 
with t' /t. Smin assumes values compatible with the ex- 
perimental findings only for (small) values of t' /t such 
that the Fermi surface is closed, while an open Fermi sur- 
face pushes Smin to exceedingly large values. This anal- 
ysis suggests that the (1^2 _y2 + id^y pairing may likely 
occur in LSCO, where the Fermi surface is expected to 
be closed, at least in the optimal and overdoped region. 
Other physical effects must be invoked for open Fermi 
surface materials like YBCO or Bi2212. 



B. Bilayer Splitting 

Materials like YBCO and Bi2212 have relatively large 
values t'/t > 0.25, and therefore open Fermi surfaces, 
which favor, according to the previous analysis, a pure 
dx2-y2 pairing. However, the strict two-dimensionality 
assumed up to now overestimates the effect of the VHS, 
which strongly favors the dx2-y2 wave. YBCO and 
Bi2212, as opposed to LSCO, are indeed multi-layer 
cuprates. Introducing the multi-layer structure of Bi2212 
and YBCO results in a "regularization" of the VHS. 
ARPES measurements in Bi2212 show indeed a bilayer 
splitting 2t± at the M points which, although smaller 
than the one predicted by band-structure calculations^Si, 
is sizable, and varies from 88 uieV to 140 meV— . We 
are not aware of similar data on YBCO (where ARPES 
measurements are much harder to perform), but the bare 
values of t± in Bi2212 and YBCO are similar, so that we 
can expect similar renormalized values. 

To analyse the role of a bilayer structure we introduce 



a band splitting of the form ^t^{coskx — cosfcj,)^ which 
modifies the effective DOS and splits the VHS partially 
spoiling the dx2_y2 principal component. On the other 
hand, the nodal regions are not affected and so one ex- 
pects the ratio V^^ /Vi to decrease as soon as t±_ / 0. 
In Fig|21we show the doping dependence of y^' /Vi for 
the same values of Vi and of cjq as in the case of Fig^ 
The two curves with open and solid squares corresponds 
to the case t±_ = O.li and 0.2t respectively and are, as 
expected, lower than that for t^^ — Q. As can be seen in 
the figure, the second of these curves, which still compat- 
ible with the measured value of t± is Bi2212, has been 
lowered enough to cross the V^"^ /Vi — 1 line. 

In conclusion, the bilayer splitting helps the dxy com- 
ponent to establish at least in all those doping regions 
for which the chemical potential falls in the dip of the 
density of states centered at — 4i', leaving one of the two 
Fermi surfaces open. This is a relevant scenario for the 
case of YBCO and Bi2212 in which the Fermi surfaces is 
open. 



C. Orthorhombic Distortion 

Another important effect that lowers the impact of the 
VHS on the pairing in the cuprates is an orthorhombic 
splitting of the bands. This is accounted for by taking 
a dispersion of the form = ^2(ta; cos -t- ty cos ky) + 
At' cos kx cos ky — n and tx 7^ ty. In this way and the 
kernel containing iSk = [^k + l^kP]^^^ that enters the 
equations (12)),®, are no longer symmetric functions with 
respect to the group. In any case this kernel can still 
be expanded in terms of the C4„-harmonics and one can 
solve the resulting coupled equations. However, if the or- 
thorhombic distortion is small, one can neglect the effects 
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on the symmetry of the gap and take only the larger ef- 
fects on the density of states. In doing this, one neglects 
the (i+-s+ and d^-s^ mixing that arises because the C4„- 
group symmetry of the system has been lowered to C2v 
In fact, the basis functions for the two one-dimensional 
representations of C2v that couples to the singlet are given 
by '''^r k~'~^r k respectively, and what we do in the follow- 
ing is to project our solutions on the original d^^-y^ and 
dxy harmonics in order to understand the role of a small 
orthorhombic distortion on the d^^^y-i + id^y solution. 

The density of states is then qualitatively similar to 
the bilayered one. In fact the VHS is split and one ex- 
pects that when the chemical potential lies in the dip the 
secondary dxy component is again favored with respect 
to the case with tx = ty. The results are represented by 
the two curves with crosses and stars in Fig|21 We have 
used for tx the same value of 200 meV as before for t and 
varied the ratio {ty — tx)/tx denoted as 7. Vi and loq are 
still those of the FigQ] In the case of 7 = 0.1 (crosses) 
the ratio V^"^ /Vi, although lower than for the isotropic 
one, is still greater than 1, as well as that for 7 = 0.2 
(stars) even if in this last one the characteristic upward 
curvature coming from the presence of the dip in the DOS 
can be clearly seen. As for the bilayer case we have used 
a constant value for Vi aiming at analysing the effects of 
band structure of the orthorhombic distorted compounds 
on the critical coupling, and we have evidenced that the 
dxy secondary component is favored by this effect. The 
values of 7 that we have used are representative of the 
case of the cuprates according to, e.g., Refi^ where a 
value of 7 ^ 0.1 is found by means of ab-initio calcula- 
tions for the case of LSCO. 



D. Quantum Criticality around Optimal doping 

We now discuss how the form of the pair potential can 
be important to obtain the dx^ _y2 +idxy order parameter. 
In particular, we discuss how this phenomenology can be 
related to the QCP scenarios for HTCS22i^2i2i. In this 
context, it has been suggested that the pairing may be 
mediated by quasi-critical fluctuations close to the QCP 
in the charge and/or spin sector(s) located near optimal 
doping. In this framework the correlation function of the 
critical fluctuations determines the functional form of the 
pairing potential. A few different QCP have been pro- 
posed, ranging from magnetic and/or charge ordering, to 
time reversal breaking. In the last part of this section we 
briefly discuss the case of charge and spin ordering, while 
here we first discuss some general consequences of a QCP 
near optimal doping which are independent on the nature 
of the ordered phase. Far from criticality (i.e., far from 
optimal doping) the short-range components Vi,V2, ■ ■ ■ 
may decrease strongly with the distance. At critical- 
ity the correlation saturates below a non-universal length 
scale, and all the components of the potential below this 
scale become of the same order of magnitude. V2 then 
becomes of the order of Vi, favoring the appearance of 



a secondary component of the gap. This establishes the 
connection between the dx2-y2 — > dx2_y2+idxy transition 
and the underlying quantum criticality. One may notice 
that, close to the QCP, where various couplings are of the 
same order of magnitude, higher-order harmonics should 
in principle be considered within each irreducible repre- 
sentation. The main outcome of this extension would 
only be a more detailed description of the gap profile, 
without changing the symmetry of the order parameter. 
Observed deviations from a pure dx2-y2 gap can be ex- 
plained in this context^^. 

We now consider the specific case of a pairing interac- 
tion mediated by the fluctuations around a charge or- 
dering QCP, with an incommensurate ordering vector 
Q = [QxiQy)- The effective interaction between quasi- 
particles in the Cooper channel can be written as 

where J7 is a residual local repulsion, Vc is the attrac- 
tion strength, is a "mass term" measuring the dis- 
tance from criticality, and Wkk'(Q") — 2(2 — cos{kx — 
k'x — Qx) — cos{ky — k'y — Qy)) contains the momentum 
dependence of the interaction. The sum over the eight 
momenta Q" = {±Qx, ±Qjf), {^Qy, ^Qx) makes the in- 
teraction symmetric under the C41, group. 

The real-space expression for (0J) is particularly useful 
in light of the previous analysis. We can in fact write 

y(r) = USr,o - ^V,ariK^)w'^^. (5) 

The local repulsion obviously enter only in the on-site 
coupling and does not therefore play any role in the 
dx2^y2 — > dx2^y2 + idxy transition. UriK^) are dimen- 
sionless functions of and the lattice distances r. When 
the system is at the critical point, i.e., in the limit 
^ 0, all the ckr approach the same value regardless 
of the index r, as we expected from the general argu- 
ments given above. The dependence on the critical vec- 
tor Q is entirely contained in the symmetry factor w^.q, 
i.e., in the above defined harmonics computed at the or- 
dering vector Q. The first two coupling constants that 
enter the BCS equations l(2Jl and Q are then given by 
Vi = ^Kaiiy^.q = ^Vcai{cosQx + cos Qy) and V2 = 

^Vca2W2.Q — Vca2C0sQxC0sQ^. The modulation of 
the couplings introduced by the critical momentum Q 
determines regions in which some coupling may become 
repulsive, thus making the correspondent harmonic im- 
possible to establish. In particular, the dj.2_y2 -|- idxy 
symmetry is possible only if both Vi and V2 are attrac- 
tive, i.e., for Qx G [-7r/2 : 7r/2] and Qy € [-7r/2 : n/2]. 
Moreover, the ratio V2/V1 is proportional to the form 
factor 2 cos Qx cos Qy / {cos Qx + cos Qy), and it is there- 
fore maximum at Q = (0, 0). The mixed symmetry order 
parameter is therefore more likely for small ordering vec- 
tors. 
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FIG. 3: Phase diagram for d^2_y2 and d^2_y2 + id^y pairing 
in the HTCS. The parameters, which are appropriate for open 
Fermi surface cuprates, are given in the text. 



The contribution to the pairing interaction coming 
from the spin sector can be written in the same way 
as with an opposite sign. This implies that 

pairing mechanisms based on antiferromagnetic spin 
fluctuations^, with ordering vector Q = (7r/2,7r/2) 
strongly favor dx2_y2 pairing, which is maximum for this 
vector, while the V2 coupling is repulsive for the same 
vector. Even taking into account small incommensura- 
tions around the antiferromagnetic vector, spin fluctu- 
ations mechanisms do not favor the appearance of the 
da;2_j,2 + idxy order parameter. Because of the relevance 
of antiferromagnetic spin fluctuations in the underdoped 
region of the phase diagram of the cuprates this analysis 
indicates that the mixed d^^-y^ -\-idxy state should not be 
present at least in the extremely underdoped compounds, 
in agreement with the experiments. 

To summarize the main outcome of this section, the 
proximity to a QCP helps in meeting the necessary con- 
ditions for the mixed symmetry (i.e., strong coupling and 
V2/V1 close to 1). Nevertheless, the specific nature of the 
QCP (like, e.g., charge and/or spin ordering) may intro- 
duce sign modulations of the real-space interaction, mak- 
ing some components repulsive, and therefore suppress- 
ing the correspondent harmonics. The pairing mech- 
anism based on antiferromagnetic spin fluctuations'^^ 
strongly favors the dx2_y2 pairing, but suppresses the dxy 
channel, leading to pure dx^_y2. Mechanisms based on 
charge fluctuations support instead a secondary dxy com- 
ponent besides a primary dx2^y2 provided the character- 
istic ordering wave vector Q is not too large (|Q| < tt/2). 

III. PHASE DIAGRAM 

We can now draw a phase diagram for the dx^~y^ +idxy 
pairing by fixing Vi for each doping S, such that the 



mean field critical temperature for dx^-y^ pairing (which 
we interpret as the temperature for Cooper pair forma- 
tion without phase coherence in the underdoped regime) 
grossly follows the behavior of the pseudogap temper- 
ature T* for small values of 6, and becomes negligibly 
small for the doping at which Tc vanishes in the over- 
doped regime. In this way we phenomenologically in- 
troduce the doping dependence of the principal compo- 
nent of the interaction Vi{S). In particular Vi{S) will be 
larger at low doping and decrease toward weak coupling 
in the overdoped regime. Around optimal doping we ob- 
tain Vi/t ~ 2, which is larger than the value we used 
to draw the diagrams of Figs. 1 and 2, and, according 
to the previous discussion on the relevance of the strong 
coupling, is more favorable for the dx^-y^ + idxy pairing. 

On the other hand, we consider the less favorable situ- 
ation in terms of bandstructure by using the single-layer 
isotropic dispersion {t± = and 7 = 0). According to the 
previous analysis, in order to have a chance to observe a 
dxy component in this unfavorable situation, we need to 
assume a small value of t' = 0.0625t which gives a closed 
Fermi surface in the relevant doping range. The value of 
the interaction is obviously crucial. Here we take, at any 
doping, V2/V1 = 0.8, a relatively large value, which, as we 
have seen, can be reasonably reached if a QCP is located 
close to optimal doping. This constant value of V2/V1 
allows us to analyse the combined role of the density 
of states and of the doping dependence of the principal 
coupling. In Fig. |31 we show Tc for dx^-y^ pairing (top 
curve) described above, and the critical temperature for 
the onset of the secondary pairing (bottom curve), cal- 
culated by solving the coupled BCS equations 10) and 
As shown in Fig. |31 the mixed dx^-y2 + idxy order 
parameter is stabilized in a small dome slightly shifted in 
the overdoped region, 0.19 < 6 < 0.24, and the dxy gap 
is always only a small fraction of the dx^-y^ gap. 

The shape of the mixed state region is the consequence 
of the balance between two main effects: i) the dxy can 
establish only when the VHS is far enough from the chem- 
ical potential, and therefore appears only above a given 
value of the doping, ii) Further overdoping, Vi{S) de- 
creases and pushes the system toward weak coupling, 
where the mixed order parameter ceases to be stable for 
reasonable values of V2/V1. The dome shape is therefore 
an unavoidable consequence of taking a sensible value for 
the ratio V2/V1 < 1. Only by releasing this physical con- 
dition, and taking V2/V1 ^ 1 the dxy component may 
indefinitely increase by overdoping, eventually leading to 
a pure dxy pairing in the heavily overdoped regioniS. The 
re-analysis of the LSCO data^^ with a dx2-y2 +idxy gap 
parameter, clearly displays the closure of the dome in 
the overdoped region, in qualitative agreement with our 
analysis. 

We note that small variations of the ratio V2/V1 can 
significantly modify the amplitude of the secondary com- 
ponent. At doping 6 = 0.22, for example, by slightly 
enhancing V2/V1 from 0.8 to 0.9, the critical tempera- 
ture for the dx2^y2 -t- idxy state is enhanced by a factor 
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of 3. 

Up to now, we have found that the dri.2_y2 —* d^2_y2 + 
idxy transition can be found for values of the parame- 
ters appropriate to the single layer isotropic case with 
open Fermi surface. The discussions of Sec. II B, C 
have shown that, as soon as small perturbations of the 
pure two-dimensional bandstructure are considered, the 
constraints for the onset of the secondary component are 
partially released. The data of Fig. 2 show in fact that for 
a realistic value of the bilayer splitting, the whole region 
0.14 < S < 0.26 becomes compatible with a dx2_y2 +idxy 
order parameter. Analogously, a small orthorhombic dis- 
tortion lowers the values of the critical coupling. Once re- 
alistic bandstructure parameters involving the above ef- 
fects are used, the appearance of the time-reversal break- 
ing order parameter in LSCO and YBCO, as well as the 
doping dependence, can be easily explained in the sim- 
ple BCS approach. The physical mechanisms leading to 
the mixed order parameter here presented have however 
a general content which does not depend explicitly on 
the use of the BCS approach. It is important to stress 
that the 'dome behavior' of the mixed order parameter 
does not depend on the details of the calculations, and 
represents a generic result within our approach. 

The experimental studies which motivated our work 
identified the transition to the mixed dx2_y2 +idxy wave 
with the splitting on the zero bias conductance peak. 
Here we finally discuss the possibility to provide a fur- 
ther experimental signature of the mixed order parameter 
by measuring the variation of the specific heat coefhcient 
lev — cv/T at dx2^y2 —> d,j,2^y2 + idxy transition. For 
the parameters used to draw the phase diagram of Fig. 
3, we have computed the mean-field jump of the specific 
heat a,i 5 — 0.22, where the secondary component is max- 
imum, and the effect of the second transition should be 
more evident. We actually computed the entropy as a 
function of temperature, and extracted the specific heat 
from cy = TdS/dT. 

In this specific case, at the transition from the metal to 
the dx^-y^ superconductor we obtain A'jcv/lcv = 1.00. 
At the dx2-y2 +idxy transition (which occurs approx. at 
10 K in this case) we obtain instead A'jcv/lcv — 0.28. 
We note that the absolute value of A^cv in the latter 
case is about two orders of magnitude smaller then in 
the former, higher temperature transition. Even though 
fluctuation effects are expected to round the mean-field 
jump, our estimated value of A^cv suggests that the ef- 
fect could be still experimentally detected. 

IV. CONCLUSIONS 

We have performed a systematic analysis of the possi- 
bility of a dx^-y2 + idxy symmetry pairing in the cuprates 



by means of the simple BCS approach which is indeed 
justified by the experimental evidence that the mixed or- 
der parameter establishes in the overdoped region. Our 
results are compatible with all the available indications 
of d^2_„2 + id xy pairing in YBCO and LSCO. We have 
separately discussed the role of the bandstructure effects 
and of the pairing potential. Since the VHS is an ob- 
stacle to the mixed pairing, all the effects that lower the 
impact of the singularity result in an enhanced tendency 
toward the dx^-y^ + idxy pairing. In particular, the bi- 
layer splitting characteristic of multilayer cuprates and 
the orthorhombic splitting of the bands both favor the 
mixed order parameter. The relevance of the bilayer 
splitting suggests that, Bi2212 is also likely to show a 
mixed order parameter, with a small mixing of the same 
order of YBCO. 

As far as the pairing potential is concerned, we have 
shown that a mixed-symmetry superconducting state is 
more likely to establish at strong coupling. In particular, 
it can appear for V2/V1 < 1 only at strong couphng. The 
presence of a QCP around optimal doping is an impor- 
tant ingredient, since it enhances the coupling and makes 
the next-nearest neighbor component of the pairing po- 
tential (responsible for the dxy pairing) of the same order 
of the nearest neighbor one (responsible for dx2-y2). 

The main result of our analysis is that a mixed order 
parameter of the form dx2^y2 + idxy can be stabilized 
in a small dome contained in the larger dx2^y2 region 
with a pairing potential being a decreasing function of the 
distance, i.e., V2 < Vi. Contrary to other scenarios, the 
dxy component does not grow indefinitely with increasing 
doping and it is always smaller than the dx^-y^. The only 
way to obtain an indefinitely increasing dxy component 
is to take an unphysically large ratio V2/V1 » 1. 

We have also shown that dx2-y2 + idxy pairing is more 
likely to occur for materials with closed Fermi surfaces 
such as LSCO, and in compounds where an inter-layer 
hopping splits the VHS. The less important role of the 
VHS in the electron-doped materials may lead to the 
dx2^y2 +idxy state also in these family of compounds. We 
have also discussed the impact of the secondary transition 
on specific heat measurements, showing that the relative 
jump of the linear specific heat coefficient is a fraction of 
the jump at the primary (higher temperature) transition 
to the dx2-y2 superconducting phase. 
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